APS/123-QED 


Vacancy Mediated Mechanism of Nitrogen Substitution in Carbon Nanotubes 

Deepak Srivastava* 

NASA Ames Research Center , CSC , Mail Stop T27-A1, Moffett Field, CA 94035-1000 

Madhu Menon^ 

Department of Physics and Astronomy , University of Kentucky, Lexington,- KY 40506-0055 
Center for Computational Sciences, University of Kentucky, Lexington, KY 40506-0045 

Bindu Sadanadan and Apparao M. Rao* 

Department of Physics and Astronomy , Clemson University, Clemson, SC 29634 

(Dated: May 23, 2003) 

Nitrogen substitution reaction in a graphene sheet and carbon nanotubes of different diameter 
are investigated using the generalized tight-binding molecular dynamics method. The formation of 
a vacancy in curved graphene sheet or a carbon nanotube is found to cause a curvature dependent 
local reconstruction of the surface. Our simulations and analysis show that vacancy mediated N 
substitution (rather than N chemisorption) is favored on the surface of nanotubes with diameter 


larger than 8 nm. This predicted value of the 

confirmed by experimental results presented. 

PACS numbers: 61.48.+c,71.15.Pd,81.10.Bk 

The discovery of carbon nanotubes (CNTs) by 
Iijima,[l] has set off intensive studies of these quasi one- 
dimensional structures. The single- wall nanotubes con- 
sist of rolled-up single graphene sheets of carbon atoms 
into nanotubes of different diameters and chiralities. Due 
to their chirality dependent electrical properties, and ex- 
cellent mechanical and thermal properties, CNTs have 
been suggested for applications in molecular electronics, 
functional composite materials, thermal management on 
a chip, and gas or drug storage and delivery devices. [2, 3] 
The possibility of the existence of nanotubes made of 
non-carbon elements such as boron-nitride(BN), B and 
N doped carbon nanotubes have also received signifi- 
cant attention, [4-1 6] where successful synthesis involve 
arc discharge methods, [8], reduction of carbon nanotube 
via substitution reactions, [10-1 2] and thermal CVD,[16] 

Significance of substitution reaction in a carbon nan- 
otube by B, N and other species is considerable from both 
materials and electronics perspective. Multiwall N doped 
CNTs have been synthesized in large quantities with the 
morphology of a bamboo-shoot like structure, [16] and B 
and or N doping of C nanotube can be used to engineer 
the band gaps of these quasi one dimensional materi- 
als for electronics applications. [17] Moreover, specifically 
substitutional doping of carbon nanotubes with B and N 
type elements has also been reported at high-temperature 
annealing experiments. [11, 12] By definition, a substitu- 
tion reaction of a N atom in a CNT involves replacement 
of C atom in the structure of the CNT with an adsorbing 
N atom in the surroundings. Therefore, conceptually a 
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critical minimum diameter for N incorporation is 


mechanism for substitution reaction with in a nanotube 
may involve a three step process: (a) creation of a va- 
cancy or a defect in an otherwise perfect structure of the 
CNT, and (b) incorporation of the external substituting 
atom (like N in this case) at the newly created vacancy 
or defect site, and (c) equilibration or redistribution, via 
diffusion, of one or many substituted atoms into energeti- 
cally favorable morphologies or configurations in the sys- 
tems. Activation barriers involved in any or all of these 
three contributing steps may be high but much of the 
synthesis processes are performed at high temperatures 
and extreme conditions in which diffusion or creation of 
vacancies may be very facile. 

In this Letter, we investigate vacancy mediated mech- 
anism of N doping of C nanotubes through, (a) energet- 
ics of vacancy formation in the CNT wall as a function 
of nanotube diameter or curvature of the nanotube sur- 
face, and (b) incorporation of N at the vacancy sites, 
through a generalized tight-binding molecular dynam- 
ics (GTBMD) scheme. [18] The theoretical calculations 
are complemented by experimental measurements for the 
same systems. For small diameter (or high-curvature) 
CNT surfaces, the formation of a vacancy causes a re- 
construction at the vacancy site, whereas in large diame- 
ter nanotubes with low surface curvature dangling bonds 
in a pyridine like configuration are stable. We find that 
N incorporation as a function of surface curvature is of 
chemisorption type in the former and of substitution type 
in the latter. A critical lower limit of the diameter of 
CNTs for N substitution is predicted and verified by ex- 
perimental observation of diameter distribution of the N 
doped CNTs synthesized using thermal CVD. 

The theoretical methods used in the present work is 
the GTBMD scheme of Menon and Subbaswamy,[18] 
which allows for dynamic relaxation of the system with 
no symmetry constraints. Successful applications of the 
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GTBMD method to C,[18] Si[19], BN [15, 20] and SiC[21] 
systems have been reported elsewhere. The GTBMD 
method has been used to predict that zigzag BN nan- 
otubes are more stable and, therefore, should be more 
readily made in experiments. [15] Recent experimental 
findings have provided a striking confirmation of this pre- 
diction by finding that, regardless of the experimental 
conditions used in their production, zigzag arrangement 
always dominates BN nanotubes. [12] In the present work 
we use the GTBMD method to obtain structural and 
energetic characterization of vacancy formation and N 
incorporation in CNTs as a function of nanotube surface 
curvature. Our experimental method for preparing N 
doped CNTs is based on a CVD method which produces 
high quality multi-wall carbon nanotubes (MWNTs) in 
bulk quantities. [22] In this method the average tube di- 
ameter ranges between 25-30 nm, with a noticeable frac- 
tion of tubes with diameter below 10 nm.[23] 

The energetics of vacancy creation and N incorporation 
in CNTs were done by computing the total energy differ- 
ences using the GTBMD method on (i) a bare graphene 
sheet, (ii) graphene sheet with a vacancy and relaxation 
induced reconstruction and, (iii) an N atom incorporated 
at the site of the vacancy. A full set of the energetics as 
a function of curvature was then obtained by repeating 
the above sequence and curving the graphene sheet with 
curvature corresponding to the range of external diame- 
ter of N doped MWCNTs observed in experiments. This 
allows us to do many energetics calculations over a range 
of curvatures corresponding to diameters varying from 1 
nm to 80 nm to a completely flat graphene sheet. The 
sheet is chosen large enough so that the simulated region 
is sufficiently far from the edge carbon atoms. For a given 
curvature the edge C atoms are held fixed and all other 
atoms are dynamically allowed to relax to a minimum 
energy configuration. 

The strain energy as a function of surface curvature for 
a graphene sheet with and without a vacancy is shown 
in Fig.la. The energies are given relative to the cor- 
responding flat graphene sheet with and without the va- 
cancy, respectively. In the absence of a vacancy the curve 
represents the strain energy due to the surface curvature 
alone (also a measure of nanotube formation energy in a 
different context) and increases as the surface curvature 
increases (or as the tube diameter decreases) from that 
of a flat graphene sheet. As expected, the increase in the 
strain energy is monotonic for the defect free CNTs of 
smaller diameter. The results were computed for tube 
diameters as large as 80 nm, but only shown for up to 
40 nm because the changes in the energies for tube di- 
ameters larger than 30 nm are negligible at the scale of 
energetics shown in Fig.la. 

The change in the strain energy as a function of diam- 
eter for the graphene sheet containing vacancy, however, 
is not monotonic and shows an interesting low energy 
feature for small diameters (large curvature) CNTs as 
seen in Fig.la. Note that the two smooth curves shown 
in the figure are drawn over many computed points used 


in the interpolation. For clarity in discusssion we have 
shown the positions of only points A, B and C along 
the potential well curve for the graphene sheet contain- 
ing vacancy. These points correspond to diameters 2nm, 
5nm and 8nm, respectively. As the figure shows, the 
strain energy remains fairly constant beyond point C. 
Our simulations show the vacancy induced reconstruc- 
tions to undergo a qualitative change while crossing the 
point C corresponding to a nanotube diameter of 8 nm 
indicating the existence of a critical diameter. The recon- 
struction for diameters less than 8 nm is shown on the 
left panel of Fig. lb where the reconstruction correspond- 
ing to point B is shown. As seen in the figure, a C atom 
moves towards the location of the vacancy and forms two 
pentagons sharing a vertex. The central C atom in this 
case has 4 strained C-C bonds. The reconstruction for 
diameters greater than 8 nm is shown on the right panel 
of Fig. lb where the reconstruction corresponding to 
point C is shown. Here the C atoms surrounding the va- 
cancy do not move towards the location of the vacancy. 
All three C atoms around the location of the vacancy re- 
main equivalent and one can conceptually visualize this 
reconstruction as the pyridine like structure with three 
incomplete pentagons with possibly three strained C-C 
bonds in the lattice. This reconstruction remain stable 
for larger diameter nanotubes. The structure shown on 
the right panel of Fig. lb (for large diameter nanotubes, 
corresponding to point C) is, thus, energetically less sta- 
ble than the structure shown on the left panel of Fig. lb 
(corresponding to point B) due to different bonding re- 
construction. Furthermore, even though bonding recon- 
structions at A and B are the same, from Fig.la, it can 
be seen that the structure corresponding to point A is en- 
ergetically less stable than the structure at B due to the 
larger curvature induced strain at A. An interesting fea- 
ture of the reconstruction shown on the left panel of Fig. 
lb is that due to the pulling in of C atom at the shared 
vertex of the two pentagons, the structure is asymmetri- 
cally non-planar. The reconstruction shown on the right 
panel of Fig. lb, on the other hand, keeps the structure 
planar throughout. This suggests that pyridine like re- 
construction is stable on surfaces corresponding to small 
curvature or large diameter (> 8nm) nanotubes. 

The approach and incorporation of a N atom on the 
two types of reconstructions are shown in Figs. 2a and b. 

A number of trajectories were used in our simulations 
of N atom approaching the two types of reconstructions. 
We find that for off-center approaches, N atom on the 
reconstruction shown on the left panel of Fig. lb has a 
tendency to form a chemisorbed type bond with the lat- 
tice. On the reconstruction shown on the right panel of 
Fig. lb, in all cases, the N atom gets incorporated in the 
lattice as in a substitution reaction. The surface curva- 
ture and off-planar nature of the reaction in the first case 
is visible in Fig. 2a, while in Fig.2b the structure before 
and after the reaction and incorporation remains planar. 
The above simulations and analysis show that vacancy 
mediated N incorporation with substitution type mor- 
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phology will be favored for surfaces of CNTs of diameter 
of about 8 run or larger. For smaller diameters CNTs, 
initial chemisorption is favored which could be followed 
by, (a) diffusion of the N atom on the the CNT surface, 
(b) reverse reaction and also sometimes, (c) incorpora- 
tion of the N atom in the lattice with the substitution 
morphology. 

The method for producing bulk quantities of high- 
purity aligned MWCNT films through the catalytic de- 
composition of ferrocene-xylene mixture at a tempera- 
ture of 700°C has been discussed in detail elsewhere. [22] 
Briefly, a feed solution of ferrocene-xylene (0.75-1.0 at.% 
Fe) is injected into a two-stage tubular quartz reactor us- 
ing a syringe pump. High-purity aligned MWCNT films 
were produced in bulk quantities on bare quartz sub- 
strates placed inside the reactor. For producing films of 
nitrogen-doped nanotubes, melamine was introduced si- 
multaneously (C 3 H 6 N 6 , Aldrich Chemicals, T su bu ma tion 
= 300°C) during the growth of MWCNT films. A quartz 
boat containing 3 grams of melamine was placed inside 
the quartz tube reactor in a region between the pre- 
heater and the furnace. [22] The location of the quartz 
boat was carefully chosen to ensure no sublimation of 
melamine as the pre-heater and furnace attain their oper- 
ating temperatures of 200 and 700° C, respectively. After 
the pre-heater and the furnace temperatures were sta- 
bilized, the syringe pump was turned on to inject the 
ferrocene-xylene mixtures, and the entire quartz tube was 
pushed towards the exhaust side so that the boat now re- 
sides in a region of the furnace where the temperature is 
sufficiently high to sublime melamine. The ferrocene- 
xylene injection rate was maintained at 1 ml/hr, and 
upon entering the pre-heater, the ferrocene-xylene mix- 
ture was immediately volatized and swept into the fur- 
nace by a flow of argon (675 seem) and hydrogen (75 
seem). Typical run time was 15 minutes during which 
time all melamine sublimed. After the reaction, the pre- 
heater and the furnace were allowed to cool to room 
temperature in flowing argon. Carbon deposits formed 
on the walls of the quartz tube and substrates were ex- 


tracted and characterized using transmission electron mi- 
croscopy (TEM). For measuring the diameter of N-doped 
nanotubes, a cluster of nanotube deposit was dispersed 
ultrasonically in tetrahydrofuran (THF). A drop of the 
suspension was placed on a TEM grid and subsequently 
analyzed. From the TEM images the diameter of the 
nanotubes was measured and number of nanotubes as a 
function of diameter is shown in Fig. 3. From the plot it 
is evident that majority of N-doped nanotubes (> 70%) 
are in the range of 20-50 nm in diameter. Significantly, 
we do not observe any N doped CNT in the sample where 
the nanotube diameter is about 8 nm or smaller. This is 
in good agreement with the above described theoretical 
prediction where vacancy mediated substitution reaction 
of dopant species through pyridine like reconstruction is 
not favored for smaller diameter single or multiwall car- 
bon nanotubes. It is interesting to note that several re- 
ports on B -substitution or B combined with N substi- 
tution in SWNTs can be found in the literature while 
no hard evidence for N substitution in SWNTs has been 
reported. Furthermore, consistent with findings of this 
study, several studies on N substitution in MWNTs have 
been reported. [16, 17] 

In summary, a vacancy mediated substitution of N in 
carbon nanotubes as a function of tube diameter has been 
investigated and described in this work. We find that, as 
a function of tube diameter, the formation of a pyridine 
like vacancy defect and the substitution of the dopant 
species (N in our case) at the defect site is favored over 
chemisorption for tubes of diameters larger than about 
8 nm. In the smaller diameter nanotubes (< 8 nm), 
chemisorption or formation of fused pentagon type de- 
fects is favored. 
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FIG. 1: (a) Strain energy of a curved graphene sheet with and without a vacancy as a function of curvature expressed as 
corresponding nanotube diameter in nm. (b) The left and right panels show the surface reconstructions formed at the vacancy 
site for nanotube diameters corresponding to points B and C, respectively. The edge C atoms held fixed during the simulations 
are shown in red. 







FIG. 2: (a) Chemisorption type N incorporation at reconstruction site corresponding to point B in Fig. la, and (b) N 
substitution at reconstruction site corresponding to point C. 




FIG. 3: Diameter distribution of N-doped MWCNTs produced by thermal CVD experiment described in the text. 


